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Abstract 

The geometric approach to mechanics based on the Jacobi metric al- 
lows to easily construct natural mechanical systems which are integrable 
(actually separable) at a fixed value of the energy. The aim of the present 
paper is to investigate the dynamics of a simple prototype system outside 
the zero-energy hypersurface. We find that the general situation is that in 
which integrability is not preserved at arbitrary values of the energy. The 
structure of the Hamiltonian in the separating coordinates at zero energy 
allows a perturbation treatment of this system at energies slightly different 
from zero, by which we obtain an analytical proof of non-integrability. 
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1 Introduction 



The property of a Hamiltonian system to be integrable may happen to be sat- 
isfied only at a fixed value of the energy. In this case we can speak of weak 
integrability and refer to the phase-space functions which are conserved only in 
correspondence of those given energy values as weak invariants. From this point 
of view, standard integrability, with invariants which are conserved functions at 
arbitrary energies, can be referred to as strong integrability. 

Weak invariants, also called configurational invariants, have been discussed 
by Hall (1983) and by Sarlet, Leach and Cantrijn (1985). Hietarinta (1987), in 
his account of the direct methods for the search of the second invariant, also 
provides a review of all the known 2-dimensional systems admitting one or more 
configurational invariants. In the present paper we use the approach of Rosquist 
and Pucacco (1995), where quadratic invariants at arbitrary and fixed energy for 
2-dimensional systems were treated in a unified way. As shown in that paper, 
the integrability condition for quadratic invariants, corresponding to second- 
rank Killing tensors of the conformal Riemannian geometry with Jacobi metric, 
involves an arbitrary analytic function S{z). For invariants at arbitrary energy, 
the function S{z) is a second degree polynomial with real second derivative and 
the integrability condition then reduces to the classical Darboux's condition 
for quadratic invariants at arbitrary energy (Darboux, 1901; Whittaker, 1937). 
The possibility of searching for linear and quadratic invariants at fixed energy 
was also addressed and some examples of systems admitting a second quadratic 
invariant at zero energy were provided. 

Generalizing the approach, Karlovini and Rosquist (2000) have discussed the 
existence of invariants cubic in the momenta at both fixed and arbitrary energy. 
Besides giving a list of all known systems admitting a cubic strong invariant, 
they find a superintegrable system admitting a cubic configurational invariant 
related to an energy dependent linear invariant. In Karlovini, Pucacco, Rosquist 
and Samuelson (2002) we discuss the case of quartic invariants associated with 
the existence of fourth-rank Killing tensors. The results obtained in this geo- 
metric framework are in agreement with the direct method when the class of 
potentials examined is the same, as can be seen comparing the cases listed in 
the above papers with those appearing in Nakagawa and Yoshida (2001). 

The set of weakly integrable systems is very large but, on the other hand, 
they are probably of limited usefulness in physical applications. One may won- 
der however if the knowledge of a weak invariant can provide information about 
the global dynamical behaviour of the system. One possibility is to investi- 
gate the existence of integrable systems with a higher-order strong invariant 
related to the weak one (e.g. linear in the momenta, as in one example given 
in Karlovini and Rosquist, 2000). The aim of the present paper is instead 
to explore the phase-space structure of the system at energies different from 
that assuring weak integrability in the general situation in which the system is 
non-integrable. In one class of systems, we show how the property of weak inte- 
grability leads to a better understanding of the dynamics of generic systems. In 
particular, the setting in which the geometric approach casts the problem allows 
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to apply powerful tools of analytical mechanics to prove the non-integrability of 
the dynamical system. 

The plan of the paper is as follows: in section 2 we recall conditions for 
weak and strong integrability of a two-dimensional natural Hamiltonian system 
limiting the application to the existence of a second invariant quadratic in the 
momenta; in section 3 we present a class of weakly integrable systems which have 
a quasi-regular behavior at near-zero energy and for which we give a Poincare 
type argument of non integrability out of the zero energy surface; section 4 
contains the conclusions. 



2 Integrability at fixed energy 

The technique devised and applied in Rosquist and Pucacco (1995, RP here and 
in the following) , allows to find the conditions such that a 2-dimensional natural 
Hamiltonian system admits a second invariant quadratic in the momenta. In 
this section we briefly recall these results that are at the basis of the applications 
of the following sections. 

2.1 Weak quadratic invariants 

We are interested in the classical 2-dimensional systems with Hamiltonian func- 
tion 

'H=\{pI+pI) + V{x,v). (1) 

Since the Hamiltonian is time- independent, energy is conserved and motion 
takes place on the hypersurface 

H{px,Py,x,y) E. (2) 

The approach of RP amounts to find a conformal transformation to new coor- 
dinates X, Y defined by 

z ^ F{w), z = x + iy, w ^ X + iY, (3) 

generated by an arbitrary analytic function S{z) via the relation 

F'iwiz)) ^ ^/SiD- (4) 

The transformation naturally induces a canonical point transformation that 
gives the new Hamiltonian 

where it appears the conformal factor 

\S{X, Y)\ = ^Js{w)S{w) = F'{w)F'{w). (6) 
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One can then show that, if the function 

G = E-V, (7) 

the so-called "Jacobi" potential, can be expressed, in the new coordinates X, Y, 
in the form 

where A and B are arbitrary functions of their arguments, the quadratic function 

Ioipx,PY,X, Y)^\{p\^pI)- A{X) + B{Y), (9) 

commutes with the Hamiltonian Q on the surface 

n{px,Py,x,y) =0. (10) 

We therefore speak of a weak invariant of the standard Hamiltonian system, in 
contraposition with the usual notion of strong invariant by which we mean a 
phase-space function which is conserved at arbitrary values of the energy. 

The origin of the phenomenon of weak integrability in the present context 
is due to the geometrization of the dynamics via the Jacobi-Maupertuis varia- 
tional principle. For a detailed account of the Jacobi geometrization we refer 
to standard textbooks (Abraham and Marsden, 1978; Arnold, 1978; Lanczos, 
1986) and to RP where, in particular, the link is deepened between quadratic 
invariants and second-rank Killing tensors which is at the basis of this approach. 
Here we briefly recall that to the Jacobi potential (0 pertains a family of confor- 
mal Riemannian metrics parametrized by the value of the energy. Each energy 
specifies a geodesic flow on a Riemannian manifold, whose projections coincide 
with the orbits of the system defined by the standard Hamiltonian up to a time 
reparametrization. The solution of the Killing tensor equations provides an 
invariant along the given flow at a given value of the energy parameter (in par- 
ticular at zero energy). It is only the additional requirement that the solution 
of Killing equations be independent of the energy that leads to integrability at 
arbitrary energy and poses a further constraint on the generating function S{z). 
If this condition is satisfied, the function commutes with the Hamiltonian (jSJ 
on the surface (0), V-B. One can clearly also extend the investigation of higher 
rank Killing tensors or more general functional forms of the second invariant 
and to the corresponding analysis of the conditions for strong integrability. 



2.2 Strong quadratic invariants 

In RP it has been shown that, to get quadratic invariants at arbitrary energy, 
the function S{z) must satisfy the condition 

Im{S"'(z)} = 0, (11) 

that is it must be a second degree polynomial with real second derivative. Actu- 
ally, we can show that condition (|ll|l not only implies strong integrability, but 
also assures the separation of variables. We have in fact the following 
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Theorem 1: Given the conformal transformation 

dz = F'{w) dw = ^/S{z{w))dw, (12) 

condition necessary and sufficient to have 

\S{X,Y)\^As{X) + Bs{Y), (13) 

where As and Bs are real analytic functions of a single variable, is that Ullf) is 
true. 

Proof: Let us prove that implies Hll|l . Functional form is equivalent 
to say that the function \S\ must satisfy the differential equation 

\S\^XY = 0. (14) 

But 

\S\,xY = i[\S\,^^-\S\^^.^]^i[Vs{VS{y/s§),,)^^-y/^{v^{^)4ll] 

= z|5|(5,,-5,-j) =-2|5| M'S"}. (16) 

Since 15*1 is a non-null positive function everywhere, the implication is evident. 
On the contrary, Ijllfl is equal to 

S = az^ + f3z + 7, (17) 

with a real and /3, 7 complex. This implies Ijl^fl as can be verified by performing 
the coordinate transformation Q, with 

w - F-i(z) = / ^ (18) 
J ^/az^ +I3z + -f 

where H17|) has been used in defining the generating function via Q). Q.E.D. 

In RP, Sect. 4, the explicit forms of functions As and Bs are given. They can 
be of the four possible types corresponding to separability in Cartesian, polar, 
parabolic and elliptic coordinates. In this way, if we write the two arbitrary 
functions appearing in (jSJ as 

A{X-E) = EAs{X)-f{X), (19) 
B{Y;E) = EBs{Y)-g{Y), (20) 

the "true" potential takes the form 

and the Hamilton-Jacobi equation for Hamilton's characteristic function, yV(X, Y), 
associated to Hamiltonian ((Sj, takes on the explicitly separated form 

i [{W.xf + (W,y)2] + f{X) + g{Y) - E[As(X) + BsiY)] - 0, (22) 

and therefore strongly integrable two-dimensional systems with quadratic sec- 
ond invariants can be exhaustively classified. 
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2.3 Regular izat ion 

Before starting the analysis of weakly integrable system, it is useful to recall 
the fact that the conformal transformation leads, in a natural way, to a time 
reparametrization of the dynamics which can be usefully exploited in the appli- 
cations. This reparametrization has been used since a long time in the frame- 
work of the "regularization" of singular differential equations, especially in the 
applications in celestial mechanics (see, e.g. Levi-Civita, 1956; Sundman, 1912). 
We can therefore speak of a generalized regularization approach and the null 
Hamiltonian introduced below can be called on the same footing the regularized 
Hamiltonian. 

Let us suppose that we have chosen a function S{z) which generates the con- 
formal transformation which gives a system admitting a weak second invariant 
®. If we denote by the function 

Up]c+pl) + f{X)+g{Y), (23) 

which can be interpreted as the numerator appearing in the Hamiltonian |SJ) 
once the expression (|21|) is taken into account, the Poisson brackets of a generic 
phase-space function J with Tis can be written as 



{J,ns} = ^{J,Ho} - ^{J, 1^1} = U^o}-^mi^|} ^ (24) 



On the other hand, let us introduce the null Hamiltonian 

■Hu^{ns-E)\S\='Ho-E\S\. (25) 

The Poisson bracket of a function with TL^ is given by 

{J.nM}^{J.'Ha}-E{J,\S\}. (26) 

Eas. (|24(l and (|26|l are therefore equivalent with respect to conservation of J. 
The conceptual difference relates to the fact that, whereas vanishing of ea. l|26|l 
expresses the possible conservation of J in the dynamics provided by 7i_v at zero 
"energy" (i.e. TLj\f — 0), the vanishing of ca. (|^ expresses the conservation of 
J in the dynamics provided by Tis at arbitrary energy. Note that the physical 
energy E enters into Tij^j- as an arbitrary parameter. 

With the null Hamiltonian, it is automatically introduced a new (regulariz- 
ing) time variable 77 by means of 

^" = 1- ™ 

Just as the Poisson brackets H24|l express the total derivative with respect to 
the standard time 

^ = UWs}, (28) 
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the Poisson brackets H26|l express the total derivative with respect to the regu- 
larizing time 

(29) 

Correspondingly, the "regularized" equations of motion are 

^ = {x,nM}, ^^{px,n^}, (30) 

^ = {Y,H^}, ^ = {PY,H^}. (31) 
arj arj 

It is worthwhile to remark that Ti.s gives the same dynamics as the original 
physical Hamiltonian Ti.. Only different coordinates are used as a result of a 
canonical point transformation. The time reparametrization implies the relation 
between Poisson brackets associated with the original physical Hamiltonian and 
the null regularized one. Considering the phase-space function of the regularized 
system as J{px,PY, X, E) and the corresponding phase-space function of the 
original system, say I{px,Py, x, y), obtained via the recipe 

{px,PY,X,Y} {px,py,x,y}, E n{px,Py,x,y), 

the total time-derivative along the flow of these two functions are related by 

— = 151 — 
drj dt 

The conformal factor is never vanishing and therefore we extend considerations 
above on conserved quantities to functions in the original coordinate frame. 

It should also be recalled that the combination of the conformal transfor- 
mation with the introduction of the new time variable by means of 127|l can 
also be seen as a canonical transformation on the extended phase space where 
Ti. and t are new canonical coordinates (Tsiganov, 2000). 

2.4 Dynamics at zero energy 

The dynamics provided by the null or regularized Hamiltonian H25|l is particu- 
larly simple in the separable case: the two motions in the separating coordinates 
X and Y decouple and the general motion is given by a superposition with in- 
dependent arbitrary initial conditions. Clearly, the same happens in the even 
simpler case of a system constructed as above choosing an arbitrary conformal 
transformation and selecting only motions at zero energy. 

In the rj time, the equations of motion (|30|l given by Hamiltonian H23|l are 
simply 



d^X 
drj-^ 



- -nX), (32) 
= -.g'(r). (33) 
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The two "energy" equations 

^pj,+f{X) = h,, (34) 
ip|.+5(y) = h2, (35) 

with 

hi+h2 = 0, (36) 
dehmit the regions admitted to the motion by means of the inequaHties 

f{X) < /ii, (37) 
g{Y) < h2. (38) 

From the theory of conformal transformations (see, e.g., Markushevitch, 1983), 
the "isothermal" net of resulting coordinates is orthogonal (except in the finite 
number of isolated singularities where F' = 0). Since in the general case, the 
frequencies of the two motions in Ij^UI) are incommensurable, the regions ad- 
mitted by the motion, with coordinates satisfying inequalities H37|l . are densely 
filled by the representative point. This happens both in the presence of rota- 
tions and librations with the two possibilities determined by the explicit form 
of the potential functions A and B and by initial conditions. This picture is 
the same even if we describe the dynamics with the original Hamiltonian Tis- 
What changes is only the speed of the representative point along the trajectory. 

Due to the almost trivial nature of the dynamics, we would like to know if 
we can exploit this simplicity to gain information about the much more complex 
situation of arbitrary energy. In the hypothesis of integrability at zero energy, at 
energies different from zero there are two possibilities: either the integrability 
of the system is preserved or it is broken. In the first case, a second strong 
invariant exists but there is no general procedure to find it. Isolated systems 
with second invariants which are polynomial in the momenta can be identified 
with the techniques described in the papers cited in the Introduction. However, 
the general case we expect is that in which integrability breaks down. This is 
the subject of the rest of the paper, where we address the following question: 
can we use the setting developed so far to investigate and possibly predict the 
non-integrability of a given system? 

Before closing this section, we would however mention the fact that, due to 
the possibility of quite involved conformal transformations, a motion that looks 
trivial in the separating variables, may become complicated with many peculiar 
features when displayed in the physical coordinates of the original problem. 
Therefore, even the case with integrability limited to the hypersurface of zero 
energy can have direct practical applications. 

3 Systems of the class S{z) = iz^ 

We start now the analysis of 2-dimensional systems generated via general confor- 
mal transformations that do not stay in the restricted class specified by Theorem 
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1. These systems are therefore integrable (actually separable) on the zero-energy 
surface. We have selected a set of systems for which the question posed in the 
previous section can be given a meaningful answer, with the criteria of the 
simplest non-trivial situations. 
We consider the function 

Siz) = iz^. (39) 

This is the simplest polynomial which does not satisfy the constraint and 
therefore gives a potential which is not automatically integrable at arbitrary 
energy. 

Let us examine the coordinate transformation generated by H39|) . By virtue 
of (0J, the corresponding confornial transformation is given by 

F{w) = eM^w}. (40) 

Using polar coordinates so that z = x + iy = rexp(ii3) we get^ 

X = ^{§ + \nr), r=^(t?-lnr). (41) 

The X,Y ^ constant curves are respectively given by the two families 

nil}) = exp{V2X ~ ?9), r2(i9) = exp{V2Y + i?). (42) 

The range of X, Y is from —oo to +oo so that the punctured plane (with the 
origin excluded) is covered with a one-to-one correspondence. 



3.1 A class of weakly integrable systems with S{z) = iz^ 
From the treatment of section 2, it follows that the potential given by 

,w)=im^. (43) 

with the separating coordinates of (|41ll . is integrable at zero energy for arbitrary 
functions / and g. The factor r^^ comes from eq.® that, from the choice of 
eq.ljSni), gives 

|^(X,y)| =e^(^-^) =r2. (44) 

Due to the presence of the polar angle as explicit argument in combination with 
the logarithm of the radial coordinate, one must be careful in the selection of 
the arbitrary functions to avoid unphysical multivaluedness in the resulting po- 
tential function. A simple choice giving a smooth glueing of the sheets, resulting 
in a continuous and single- valued potential, is the following 

f{X) = i(C-sinV2X), (45) 
g{Y) = i(C-sinV2y), (46) 



^Note that in RP, in equations (102—104), there is a misprint and the coordinate r must 
be substituted with its natural logarithm. 
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where C is a real constant. Using (|41|l . the exphcit form of the potential in 
polar coordinates is 

Vir,^)^ ^~'^^J'<'^'\ (47) 

When the parameter C lies in the interval 

< C < 1, (48) 

the motion in the potential (|47|l is bound for every value of the energy below a 
positive threshold, since the equipotentials are closed regular curves around the 
unique absolute minimum of the function V{r, z?), 

E^,n = min{F(r, t?)} < 0. (49) 

At the threshold energy, Et > 0, the equipotential is asymptotically open and 
the motion is of parabolic type. We are interested in the bound motion, so that 
we study the energy range 

E^.rr <E<Et. (50) 

At zero energy we can readily verify that the system is indeed integrable. Using 
momenta expressed in polar coordinates according to 

rpr = xp.^+ypy, (51) 
PS = xpy - ypx, (52) 

the second invariant has the form 

Io{pr,P§,r,'d) = rprP^ — cos -ff sin(lnr). (53) 

It is straightforward to check that this function commutes with Hamiltonian 

C — sinz9cos(lnr) 



n 



1 / 2 pi\ O-smfcosllnrl 

2 [Pr + ^)+ ^2 (54) 



on the surface ifTIUl . 

Moreover, a direct inspection allows to find the main periodic orbits. We 
have the marginally stable orbit 

r = 1 (55) 
and the two stable periodic orbits given by 

i} = ±lnr. (56) 
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3.2 Dynamics outside the zero energy shell 

The interesting question is now: what is the nature of the motion when the 
energy is different from zero? Since the system has been constructed as the 
simplest generahzation of separable systems, it could be still integrable possi- 
bly with a more complicated second invariant. A systematic investigation of 
polynomial invariants of higher degree has been unsuccessful. This is clearly 
not enough to state non-integrability of system H54|l . However, we can take the 
other way and try to see if we can devise a test of non-integrability. Actually, 
it would be highly desirable to have an analytic test of non-integrability, since 
it can escape a numerical approach, such as the computation of the Poincare 
surface, if stochastic zones are exceedingly small to be detected. 

An analytic proof can be obtained by applying the Poincare theorem (Poincare, 
1892; Whittaker, 1937) on non-existence of additional invariants in non-degenerate 
systems. To this end we follow a perturbative approach in which the value of the 
energy plays the role of the perturbation parameter, with zero energy defining 
the unperturbed system. We remark that this procedure can be applied to the 
regularized system whereas it seems very difficult to attempt it in the original 
variables. 

In the present instance, the null Hamiltonian ()25|l in the separating variables 

is 

'Hj^ = Upx+Py)^ 1 (sin\/2X + sin V2r) +C-E (57) 
With the further change of variable 

a{X) = V2X- I =?9-Mnr- |, (58) 
P{Y) = V2y-|=i?-lnr-|, (59) 

we obtain the Hamiltonian 

Up = ^{pl+pl)-cosa-cosP-2Ee"-'^ = -2C. (60) 

We can now express the system in the standard form Hp = Hq + eHi, where 
an integrable "zero" order Hamiltonian is perturbed by a non-integrable first 
order term. In our case we have 

^0 = ^{p1 + pp) - cosa - cosP, (61) 
Hi = e"-'^, (62) 
e = -2E. (63) 

It is then necessary to express this Hamiltonian as a function of the action-angle 
variables of the unperturbed system: 

Hp = HoiJi, J2) + eWi(Ji, J2, 01, 02). (64) 
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The canonical variables satisfy the equations of motion 

= ^^^ajT' (65) 

where the prime denotes differentiation with respect to 77, the new time intro- 
duced in H27|l . For E = (e = 0), the unperturbed motion is given by: 

J'k = Jk ^ const, (66) 



dJk 



LUk = const 6k = LOkTi + const. (67) 



The energy equations (|34|l represent in this case the two uncoupled pendulums 
of the unperturbed Hamiltonian. Now we have 

hi + h2 = -2C. (68) 

In view of (|48|l . we are considering only libration motions. It follows that the 
action variables are (see, e.g., Boccaletti and Pucacco, 1996) 

Jfe = -[i?(4) + (4-l)X(4)], (69) 

TT 

where 

Kl = ^{l + hk), 0<4<1 (70) 

and K{Kf,) and E{k1) are the complete elliptic integrals of first and second kind. 
Ti.0 is implicitly expressed in terms of the actions by means of 

Ho(Ji, J2) = 2(k2(Ji) + k2(J2)-1). (71) 

Also the main periodic orbits of the unperturbed potential of Hamiltonian 
H6(J|I deserve to be mentioned. The two periodic orbits 

a = 0, (3 = 0, (72) 

are elliptic stable and coincide with the two periodic orbits in ea. (|56|l . The two 
periodic orbits 

a -13 = 0, (73) 

and 

a + ^ = 0, (74) 

are of parabolic type (that is they are marginally stable with characteristic 
exponents ±1). In particular H73|l coincides with the orbit of ea. H55|l . These 
orbits are "resonant" since their frequencies cui , lu2 are equal and are members 
of a family given by all the possible phase shifts A with 

< A < — . (75) 

U)l 
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Orbit H73|l corresponds to A = 0, whereas orbit H74|l corresponds to A = vr/cji. 
Parabolicity is just related to the existence of resonant tori in the unperturbed 
system. A perturbation "breaks" these tori generating isolated stable and un- 
stable periodic orbits. In general, in non-degenerate systems we may expect the 
existence of a generic resonance that in the following will be taken of the form 



with p and q integers such that p < q. 

3.3 Poincare theorem on non-existence of additional in- 
variants 

Poincare method to demonstrate non-integrability of a Hamiltonian systems is 
based on proving that, under some genericity conditions, no analytic second 
invariant exist besides the Hamiltonian itself. We state the theorem as it is 
presented in Poincare (1892) and Whittaker (1937). Then we show (theorem 5 
below) that our system H57|l comply with the requirements of Poincare theorem, 
proving therefore its non-integrability and, as a consequence, non-integrability 
of original system (|54|) . Moreover, by exploiting lemma 4 used in the proof of 
theorem 5, we will show in the following section the existence of isolated periodic 
orbits and verify that their stable (or unstable) nature varies with respect to E. 

We start with a definition. Assume that, for Ji = Ji and J2 = J2, frequencies 
(111 and UJ2 are commensurate according to relation {TSJ. Ji, J2 satisfying this 
condition are said to belong to the Poincare set. For a two dimensional system, 
"Poincare non-existence theorem" can be stated as follows: 

Theorem 2: Let the Hamiltonian 



^1 = [p/q)^2, 



(76) 



Hp — Hn{Ji, J2) + eHi(Ji, J2, 6*1, ^2) 



(77) 



satisfy the following hypotheses: 
1. Ha is non-degenerate, that is: 




(78) 



2. Hi is generic, that is, defining the Fourier expansion 




(79) 



no coefficient hmniJi, J2) is zero in the Poincare set; 
then, there is no analytic second invariant of the form 



00 




(80) 



s=0 
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independent of H. 



Proof: We refer to Poincare (1892, sects. 81-83) and Whittaker (1937, sect. 
165). For a more recent proof, see Giorgilli (2002, sect.3.1). Q.E.D. 

To show the non existence of analytic second invariants of the form H8U|I 
is equivalent to exclude Liouville integrability. However, it must be remarked 
that the above result apply in the case of a regular function which is uniformly 
continuous in the phase-space variables and in a small range of the perturbation 
parameter, but the possibility cannot be excluded that for small fixed values of 
e does a non-uniform invariant indeed exist (see, e.g., Kozlov, 1991, sect. 4.1). 

We now proceed to verify that Hamiltonian H6U|I satisfy the conditions re- 
quired to apply theorem 2. Let us check the non-degeneracy condition first. 

Lemma 3: The Hamiltonian (|71|) is non-degenerate. 

Proof: The frequencies of the unperturbed motion can be calculated as fol- 
lows (see, Boccaletti and Pucacco (1996), sect. 1.16): 

dHo d 2\ . ( dJk\^^ TT 1 



where we have used 

^-^^ \,K{4), (82) 



dKk TT 

which comes from ea. (|69f) and standard properties of the elliptic integrals. With 
an analogous procedure we can compute the second derivative 

a^Ho _ dKk d^Ho _ fdJkY^dKk _ 16 1 
84 dJkdKkdJk [dKk J dJk tt^kIK^kI)- ^ ' 

We recall that K{k'^) is a positive uniformly increasing function in the interval 
< < 1, such that 

K{0) = ^, if(l)->oo. (84) 



Since 



we have that 



det(^H^^^O, (86) 



dJl J " dJl dJl 

so that ea. (|78|l is satisfied and condition 1. in the statement of Theorem 2 is 
verified. Q.E.D. 



We come now to the genericity condition, directly proceeding to perform a 
Fourier series expansion of the perturbation appearing in H6U|I : 



+ 00 



= e'^^^^-P^^) = J2 /imn(Ji, J2)e*('"''^+"'=\ (87) 



m,n— — oo 
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where we have to use the solution of the unperturbed problem. The motion in 
the potential of Hamiltonian (|60|l with E — is given by 

a{ri) = 2 arcsin[Kisn(7y, Ki)] , (88) 
/3(r/) = 2arcsin[K2Sn(?7, K2)] , (89) 

where sn(u) is the Jacobi sine-amplitude function. In this way, since k depends 
on J, the actions of the unperturbed problem enter in the expansion. On the 
same footing, the angles appearing in the Fourier expansion are just the canoni- 
cal angle variables of the unperturbed problem, so that the harmonics mcoi, nuj2, 
with integer m, n, are multiples of the unperturbed frequencies. 

We preliminarly assess the structure of the resonating set of frequencies that 
determine the Poincare set. Using expressions (|81|) for the frequencies, we have 
that the resonance condition H76() can be written as 

-A ^ = (90) 

where, in view of and ifTUjl. 

Choosing a value of the parameter C in the interval (|48|l fixes the corresponding 
range of k\ as 

< < 1 - C. (92) 

Therefore, we can define a function 

R{kI)='^K{1-C-hI)-K{hI), (93) 
which is continuous in the range H92|l and, recalling H84|l . such that 

R{Q)^1k{1-C)~^, R{l^C)^^-K{l~C). (94) 
p 2 2p 

Without loss of generality we have assumed p < q. Recalling again H84|l . we 
have that 

which implies R(0) > 0. Moreover, if the inequality 

- < -K{1 - C) (95) 

p TT 

is satisfied, we see that R{1 — C) < 0. Then we have that the function R{k\) 
defined in the interval (|92|l . by continuity, must necessarily vanish at a point of 
this interval. Observing that condition (|9()|l can also be written in the form 

^K(1-C-kI)-K{kI) 

p— li 12: ^ f96l 
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we see that the point where R{k\) vanishes is also the point where the resonance 
condition is satisfied. 

Actually, it is straightforward to check that inequality H95|) is not a limitation 
on the set of possible resonances, since dynamics automatically comply with 
it. In fact, recalling again (|81|l and I84|) . we see that both frequencies of the 
unperturbed libration motions must satisfy the constraints 

Lo{l-C) <uj< uj{Q) = 1. (97) 



2K{1 - C) 

Therefore, the ratio of two frequencies is such that 

^<^^,=-K{l-C), VCe(0,l]. (98) 

UJl L0(1 — C) TT 

As a consequence we see that, for every pair of frequencies such that the reso- 
nance condition H76|l is verified, inequality H95|l is always satisfied. We conclude 
that, for every value of the parameter C in the range l|48|) . there is a whole 
dense set of resonances in the space of unperturbed frequencies. Each value of 
the root of the function R{k\), say k\{C,q,p) with the corresponding value of 
^2 given by l|ilT|l . determine, via (|?^ . the corresponding values of Ji and J2 of 
the Poincare set. 

At this point, we have the following: 

Lemma 4: The Fourier expansion coefficients of perturbation l|87(l are non- 
vanishing on the resonance manifold specified by relation (|76|l . 

Proof: The Fourier coefficients are given by 



where 



Am = ^ r e^^^^e-^^-i^dT?, (100) 
Jo 

i?„ = ^ r e-f^^^e-'^^^^dv, (101) 



and, in a natural way. 



Tk^—, fc = l,2. 

OJk 



On the resonance manifold we have 



TO q 
n p' 



(102) 



so that the simplest choice we can make is 

m — q, n — —p. (103) 
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Let us separate real and imaginary parts in the Fourier coefficients of a 
generic resonance. The coefficient with m = q is: 



1 /""^^ 

Aq = — e°''-^^\cos qui iT] — i sin qCji-q)dri ~ Cq — iSg. (104) 
Ti Jo 

Using the resonance condition 

the coefficient with n = —p is: 

B^p = —— / e~^^^\cosqLUir/ + ismqu}ir])dri ~ Cp + iSp. (105) 



The Fourier coefficients on the resonance manifold are therefore 

hq^^p - {CqCp + SqSp) + l{CqSp - CpSq). (106) 

We now proceed to evaluate above integrals, starting with the elimination of all 
cases in which some of them identically vanish. First of all, observing solutions 
H88H89|I . we realize that the functions a{r]) and f3{ri) are even over each half os- 
cillation period and odd over the corresponding whole period. As a consequence, 
the exponentials appearing in the perturbation (|H7jl are both even functions over 
the corresponding half period. On the same interval, cosgtDiTy (with q odd) and 
singwi?7 (with q even) are odd functions. Referring to the integrals in HlU4fl . by 
parity arguments, we can therefore state that: 

= for g odd, (107) 
Sq = Q for g even. (108) 

The remaining integrals, namely Sq with q odd and Cq with q even, are surely 
non vanishing. For, take for instance Sq when q is odd. We have 



/■--II /"Ji 

TiSq = / e"('')sinQa)ir?d7/+ / e"^'') sinqcjiryd?] (109) 

= r \°'^'^'> sinqiviTjdT]- r \~°'^'^hmqLJir]dT], (110) 
Jo Jo 



where, in the second step, we performed a translation 77 — !■ — Ti/2 and used 
the fact that 

a(7,-ri/2) = -airj) (111) 

smqui{ri — Ti/2) = —sinquiir]. (112) 

The statement follows because 
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The same holds for Cq when q is even. 

The integrals appearing in (|105|l are characterized by the integration interval 
which depends on the ratio p/q. However, the function exp(— which has 
the same parity property of exp(— a), is "sampled" by trigonometric functions 
alternatively even and odd on each half interval. Using the same arguments as 
above we get 

Cp = 0, Sp^O, for p odd, (113) 
Sp = 0, Cp^O, for p even. (114) 

Summarizing these results and comparing with expression (|106|) . we get fi- 
nally that the following coefficients 

CqCp for even, (US) 

iCqSp for q even and p odd, (116) 

—iSqCp for q odd and p even, (117) 

SqSp iovq,podd, (118) 

are non vanishing. This proves the lemma and shows that condition 2. in the 
statement of theorem 2 is indeed satisfied. Q.E.D. 

We remark that, according to the analysis already carried out by Poincare 
(1892, sect. 83), condition 2. can be somewhat relaxed. In particular, introduc- 
ing the concept of "classes" in the resonance manifold, the first combination 
with both q and p even can be reduced to one of the others. However, this fact 
does not substantially modify the conclusions above. 

With this body of results we can now enunciate our non-integrability state- 
ment applied to the original class of systems of section 3.1. 

Theorem 5: The dynamical system generated by Hamiltonian, 

1/9 ■o'\\ C — sin z9 cos(ln r) 

n = \{pl + '^) + — ^ = E, (119) 

integrable at energy E — Q, does not admit a conserved function 

/(p„p^,r,^;C), (120) 

defined in the open set — ei < £^ < and < i? < £2, with small enough ei_2, 
for every value of C in the interval (0, 1]. 

Proof: From the discussion on regularization and dynamics at zero energy 
of sects. 2.3 and 2.4, it has been shown that the dynamics generated by Hamil- 
tonian (|119|l is equivalent, a part a time reparametrization, to those generated 
by Hamiltonian (|6()(l . In particular, this applies to the conservation properties 
of phase-space functions: non-existence of a function I{pa,Pi3, a, P; E) which is 
an invariant, for any i?, of Hamiltonian (|60|l . implies non-existence of a function 



''q,-p 
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I{pr,P'd,r,'d; C) which is an invariant, for any C, of Hamiltonian pi9|l . Hamil- 
tonian H60|l satisfies the conditions for the applicability of theorem 2. For, by 
lemma 3 we have that its unperturbed part is non degenerate and from lemma 
4 we have that the perturbing function is generic in the Poincare set which is 
dense for every < C < 1. The statement follows. Q.E.D. 

We recall (see sect. 3.1) that Hamiltonian (I119|l supports bound motion in 
the energy range 

Emin < E < Et (121) 

with Emin < and Et > and is separable a,t E — 0. Theorem 5 proves 
non-integrability of these motions for small (positive or negative) values of the 
energy. However, from the remarks above, it is not possible to exclude the 
existence of isolated integrable case with non uniform invariants corresponding 
to some other fixed values of the energy. 



3.4 Existence of isolated periodic orbits 

Non-integrability of a perturbed Hamiltonian system manifests itself with the 
birth of alternately stable and unstable periodic orbits. Actually, their existence 
is an "obstruction" to the integrability and can be used as an argument to prove 
non-integrability when the conditions required by theorem 2 do not apply or are 
difficult to be checked. In our case, we can exploit the results already obtained 
to get informations about the nature of the periodic orbits ensuing from the 
perturbation. We follow Arnold, Kozlov and Neishtadt (1988, sect. 6.1) to 
show the existence of isolated periodic orbits and to determine their stability 
properties. For a recent application of this technique in celestial mechanics we 
refer to Diacu and Santoprete (2001). 

We state, referring for the proof again to Poincare (1892), the following 

Theorem 6: Suppose that the following conditions arc satisfied: 

1. TYo is non-degenerate in J^., fc = 1, 2, that is: 

2. for some A = A*, the average perturbation, defined as 

Wl(Jl, J2,A) = ;^ f ' ■Hi{Jij2,^iV+\^2V)dih (123) 



is such that 



dX ' dX 
and the perturbing function 



''^^ 0, ?S^^0; (124) 



Til = Wi(Ji, J2,(ii?7 + A,W2?y), (125) 
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where A is the phase shift of the resonating motions, is periodic in rj with period 



Then, for small e, the perturbed system has a periodic solution with period 
T^(e) such that T'r;(0) = and with characteristic exponents ±a which can be 
expressed as power series in the form 



a — a\\fk + a-it + azi^\fk + dots, 
where the coefficients are given by the relations 

uJi ai = j^yt I ^1 „ 2ujiuj2 - 



5A2 



9J2 



' dJidJ2 



(126) 



(127) 



We can actually prove that system l|f)l|l meets the requirement of theorem 
6. For, condition 1. is a straightforward application of lemma 3 to the Poincare 
set. Condition 2. derives applying the results of lemma 4. Taking into account 
the phase shift, the Fourier expansion of the perturbation can now be written 
in the form 

Hi(Ji, J2,??,A) =^i„(A)i3„e^('""l+"'^=)^ (128) 

m,n 



where 



and 



imcJi A 



Ti 



i(r)- 



(129) 



r = ry + A, 

In order to verify condition 2. we have to perform the average H123|) of TYi. Using 
the expansion H128|) . we realize that the only A dependent terms surviving the 
averaging process, are just those corresponding to resonances of the form ()76|l. 
that is 

Hi(Ji, J2,ry,A) = Vi„(A)S„e'("+"f)'^i''. (130) 



Therefore, after averaging, we get terms of the form 

Hi(Ji, J2, A) = 2Re{e"?^i^A,B_p} 



(131) 



With the same notation used in the proof of lemma 4, first and second derivative 
of Hi{Ji, 0/2, A) with respect to A are: 



dTii 
dX 



2Re < iquji 



^ = 2Re^,2^2 



C^1^1^{Cg+tSg){Cp + tSp) 



(132) 



(133) 



Therefore, condition 2. is satisfied if a A = A* exists such that the imaginary 
part of the term in the square bracket vanishes with a non vanishing real part. 
The imaginary part is 



cos q(jJiX{CqSp + GpSq) + saiqCji\{CqGp — SqSp). 



(134) 
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This expression vanishes if A takes the value 



A-^^— arctan^i^±^. (135) 



Correspondingly, we have 



2q'Cji cos qujiX*[Cq{Cp + Sp) + Sq{Cp - Sp)]. (136) 



Using the results summarized in (|1(J7II and (|113|l . if n denotes an arbitrary 
integer, we get that the critical shifts are either 

A* = ^ (137) 

qu}i 

in case of an odd-odd or even-even q,p combination, or 

= ^ (138) 

in case of an even-odd or odd-even q,p combination. 

As a simple application, let us take the simplest case, that is 

p = g = 1, a)i = 0)2, 

so that the two unperturbed motions are identical but for the phase shift. From 
P37|l we then have 

A* = ^. 

In particular, 

A* = ^ ^ = '^i'^'=i > ^' 

so that, for the zero phase shift corresponding to the unperturbed periodic orbit 
H73(l , the leading order term in the expansion (|126|l of the characteristic exponent 
is a real number. This can be seen using ea. (|127|) with the result H139|) above 
and with the round bracket that is positive in virtue of l|83|) and (|85|l . Therefore, 
the characteristic exponents are real if i? < and pure imaginary if i? > 0. The 
periodic orbit is then unstable in the former case and stable in the second. This 
change in the nature of a periodic orbit on varying the energy is another typical 
feature of non-integrable dynamics. 



4 Conclusions 

The main points which we have focused on in the present paper can be summa- 
rized as follows: 

— There exist wide classes of natural conservative mechanical systems with 
integrable (actually separable) dynamics at energy equal to zero. 
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^ In the general case, the dynamics of the system is non-integrable at en- 
ergies different from zero. 

— The technique exploited here to construct these systems, which is based 
on a generic eonformal transformation, offers the possibility of characterizing 
many features of the non-integrable dynamics. 

In particular, we have seen how the structure of the transformed Hamiltonian 
naturally allows a perturbative approach to investigate the departure from the 
integrable regime. Standard tools of analytical mechanics, like the Poincare 
methods, can be applied more easily than in the usual setting, to get analytical 
proofs of non-integrability. 

A natural continuation of this work is to explore the character of other 
systems constructed by means of other eonformal transformations or even of 
other potentials in the classes introduced above. In particular, we have limited 
the analysis to systems admitting bound motion, but in the class of systems 
studied above, it is easy to construct potentials allowing the coexistence of 
limited and unlimited motions. In such systems, dynamics are characterized by a 
transition to chaotic scattering with possible applications in celestial mechanics, 
quantum mechanics and general relativity. 
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